Cherenkov radiation from fluxon in a stack of coupled long Josephson junctions. 
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We present a systematic study of the Cherenkov radiation of Josephson plasma waves by fast moving 
fluxon in a stack of coupled long Josephson junctions for different fluxon modes. It is found that at 
some values of parameters current-voltage characteristic may exhibit a region of the back-bending 
on the fluxon step. In the opposite limit the emission of the Cherenkov radiation takes place. In 
the annular junctions of moderate length the interaction of the emitted waves with fluxon results in 
the novel resonances which emerge on the top of the fluxon step. We present more exact formulas 
which describe the position of such resonances taking into account difference between junction and 
non-linear corrections. The possibility of direct detection of the Cherenkov radiation in junctions of 
linear geometry is discussed. 
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I. INTRODUCTION 

Cherenkov radiation exists if a particle moves with 
a velocity equal to the phase velocity of the emitted 
waves. Electromagnetic waves in a long Josephson junc- 
tion (LJJ) are described by the sine-Gordon equation 
which has, in particular, soliton and plasma wave so- 
lutions. Physically, the solitons are Josephson fluxons 
(often called Josephson vortices), and plasma waves in 
most cases are just small amplitude linear electromag- 
netic waves. A soliton behaves here as a quasiparticle 
with its own characteristic mass and velocity. Therefore, 
one may consider the problem of Cherenkov radiation of 
plasma waves induced by fast moving solitons. 

The dynamics of Josephson phase in a LJJ is 
rather accurately" described by the perturbed sine- 
Gordon equationEl: 



i<f>^ - 7 



(1) 



where r.h.s. of (El) is usually referred as perturbation. The 
coordinate x is normalized to the Josephson penetration 
depth A J, the time t is normalized to the inverse plasma 
frequency 1/ujp, 7 = j /jc is the bias current density j nor- 
malized to the critical current density j^ a = 1/v^ is 
the damping coefficient, /3c is the well known McCumber- 
Stewart parameter. Generally, it is assumed that a and 
7 are small. 

If we neglect the dissipation and the bias current [r.h.s. 
of (p^)] and consider a small amplitude wave 



= A siii{kx - Lui) , y4 < 1 



(2) 



by substituting (0) into (nh we obtain the well known 
dispersion relation for small amplitude plasma waves in 
LJJ: 



uj = \/h? + 1 . 
The phase velocity of these waves is 



^ph 






(3) 
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and can take any value between 1 and c». Here the ve- 
locity u is normalized to the, so-called, Swihart velocity 
Co = WpAj. From (|j) it follows that the Swihart velocity 
is the minimum phase velocity of electromagnetic waves 
in the system. 

The fluxon moving with the velocity u in a LJJ can be 
described by the following solution of (|l|) without a and 
7 terms: 
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Taking into account the r.h.s. of (|l|) as a perturbationtl 
one may derive the force-balance dependence of the 
fluxon velocity u on the bias current 7 and the dissipation 
coefficient a: 



1 



(6) 



The allowed range of fluxon velocities lies between and 
1. Thus, the maximum fluxon velocity Umax coincides 
with the minimum possible phase velocity cq of linear 
electromagnetic waves in the system and Cherenkov ra- 
diation can not take place. This coincidence results from 
the structure of sine-Gordon equation and, in principle, 
may not exist in more complex systems such as coupled 
LJJ's, arrays of JJ's or systems with special resonant 



structureaa. Coupled systems received a lot of attention 
during the last few years due to the progress achieved 
in fabrication of Nb-Al-AlO^-Nb low-Tc tunnel junctions 
and studies of layered high-Tc|-£uperconductors that ex- 
hibit intrinsic Joscphson cffectcl. 

The behavior of Joscphson phases in a system of cou- 
pled LJJ's is described by a system of coupled sine- 
Gordon equationaj. Since the dispersion relation for lin- 
ear waves as well as the maximum velocity of a fluxon 
are influenced by mutual coupling between the junctions, 
one may find the conditions at which Cherenkov radia- 
tion appears in such a system. Recently we reported clear 
evidenced of radiation by fluxons moving in a system of 2 
coupled annular LJJ's. In this paper we present system- 
atic study of the dependence of this effect on junction 
parameters, location of the trapped fluxon and report a 
more detailed model then the one presented in Ref. |^. 

The paper is organized as follows. In section II we 
present a model which extends the model of Ref. |5| in 
several aspects. The numerical technique and the results 
of simulations are presented in section III. The theory 
from section II is further extended to describe the po- 
sition of resonanaces obtained in simulation more accu- 
rately. The question whether Cherenkov radiation can be 
observed directly in LJJ of linear geometry is discussed. 
Section IV shows experimental results. The phenomenon 
of Cherenkov radiation in A^-fold {N > 2) stack is dis- 
cussed in section V. Section VI concludes the work. 
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in such a system and linearizing sine/) '^ « cj) '^ we ob- 
tain a dispersion relation which consists of two branches 
corresponding to two different modes of small-amplitude 
electromagnetic wave propagation in the system: in- 
phase mode and outjof-phase mode. In general, this rela- 
tion is rather bulkylj but for the case D' = C = 1, which 
we will be mostly interested in, it reduces to the simpler 
form: 
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From (P) the phase velocity Uph — uj(k)/k can be calcu- 
lated as 
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The phase velocities of the in-phase mode ["-I-" sign in 
(0) and (lOh] are in the range from c+ to 00 and for the 
out-of-phase mode ["— " sign in (^) and (10)] from c. 
00. In the general case, the velocities c± are defined ai 



II. MODEL FOR TWO COUPLED JUNCTIONS 

From both practical and theoretical point of view, 
it is interesting first to consider the most simple sys- 
tem such as two coupled (generally asymmetric) LJJ's, 
which we denote as LJJ and LJJ^ from now on. They 
can be described by two coupled perturbed sine-Gordon 
equations 



7^ 



(7) 

where S (—1 < S < 0) is a dimensionless coupling con- 
stant, D' — d'^ /d'^ is the ratio of the effective magnetic 
thicknesses of LJJ^ and LJJ^, C = C^/C^ is the ratio 
of specific capacitances, J — j^ / jf is the ratio of critical 
currents, a^'^ and 7"^^-^ — j^'^ / jf'^ are the damping 
coefficients and normalized bias currents, respectively. 
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and have the same meaning as the Swihart velocity Cq 
for a single LJJ but belong to different modes of electro- 
magnetic wave propagation. In the case D' = C = 1 the 
expression ( pi] ) for Swihart velocities reduces to: 
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c± = 
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I^otc, that in the this case c± do not depend on J. 

In a single LJJ, the maximum fluxon velocity Umax co- 
incides with the Swihart velocity of plasma waves cq. But 
in general, the Swihart velocities [e.g. (pi])] have nothing 
in common with the possible range of fluxon velocities 
in the coupled system since no solution similar to (o) is 
known in general case. 



B. Fluxons 



A. Plasma waves 



For the analysis of linear modes we use Eq. (|7|) in the 
absence of perturbations. Considering the small ampli- 
tude waves 



The variety of fluxon configurations is very rich even 
for 2 coupled LJJ's. When discussing different fluxon 
configurations we will use the notation [iV|M], which 
means N fluxons located (trapped) in LJJ"^ and M flux- 
ons in LJJ^ {N, M < describe anti-fluxons). 

The [1|0] state. Let's consider now two coupled LJJ 
with J 7^ 1 and some travelling wave solution (f)^'^{x — 



ut) . Introducing a new variable ^ = x — ut and omitting 
a and 7 terms, the unperturbed Eqs. ([7]) can be rewritten 
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Let's now suppose that the solution (f) '^ moves with 
velocity m = c_ (11). Substituting u ^ C- from ( |ll| ) into 
(O) we obtain the following relation between 0"* and 0^ : 
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where 



D'C - 1 + .yiD'C - 1)2 + AD'CS"^ 
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In the simplest case D' — C — 1 and J 7^ 1 we have: 
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After establishing Eq. (Mj), let's consider [1|0] state 
with the only fluxon in LJJ , i.e., (j)'^ grows from to 27r 
as X changes fr'oni —00 to +00 and 0^ changes with x but 
</)^(±oo) =0. If K = 1 then from Eq. (|l|) it follows that 
^-^ = 0^ + 27rfc. Thus, if 0^ grows from to 27r, (j)^ also 
grows in the same way which corresponds to the state 
[1|1] and contradicts to the assumption that we consider 
the state [1|0]. Therefore in the state [1|0] with k = 1 
the fluxon can not move with the velocity c_ — it is the 
maximum asymptotic (not reachable) velocity of fluxon. 

If K > 1 and (j)"^ grows from to 27r as a: changes 
from —00 to +00, 
from —00 to +00 

approaches 7r/2 at some point a; = ccq, i.e., sin(0"^) ap- 
proaches 1. At this moment sin((^^) will approach 1/k 
(k > 1) and 4>^ will be equal to some value between 
and 7r/2 {e.g. tt/S). At the next point x = xi > xq where 
(j)"^ comes over 7r/2, sin(0"^) starts to decrease, so sin(0-^) 



^ will not make 27r-leap as x changes 
Imagine that (j) grows from and 



becomes smaller than l/n and < 7r/3, i.e. 



will 



not make a 27r-leap. This is an important point: the fact 
that 6^ does not twist means that there is no fluxon in 



LJJ^ 



I.e.. 



[1|0] state is possible at the velocity u = C- 



From mathematical point of view (for a "good" differ- 
ential equation) , if a solution exists at point w = c_ it will 
also exist in the vicinity of this point, i.e., also above c_. 
In fact, simulations show that for D' — C — 1, J — 0.5, 
S = —0.5 and a = 0.1 the state [1|0] can survive up to 
about Co at some parameters of the system. This region 
of velocities c_ < u < wlnLx is a domain were Cherenkov 
radiation in [1|0] state is excited by a moving fluxon. Of 

course WmLx is a function of such parameters as J, S, a. 



For K < 1 a fluxon in the [1|0] state can not reach 
the velocity c_. Moreover, the maximum velocity of the 
fluxon in this case is smaller than c_. 

Here we have to stress that rather special conditions 
are required to get Cherenkov emission in this system. 
The most important condition is the asymmetry of the 
two-fold stack. For the [1|0] state not only the condition 
J < 1 leads to Umax > c_ but also, e.g., 7^ < 7"* results 
in a similar phenomenon as we checked numerically. If 
one includes 7"* ^ 7^ into the coupled equations not as 
a perturbation but as initial offset of (p^ , it is possible to 
obtain an equation very similar to the above equations 
for J ^1. 

The [1|±1] state. If LJJ's are identical {C = D' = J = 
a^/a^ = 1), the coupled sine-Gordon system consists of 
two identical equations and, due to the symmetry, in the 
state [1| ± 1] the relation 0"* = ±0^ can be satisfied. 
Therefore, two Eqs. (H) reduce to one equation 
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which is the usual sine-Gordon equation with the Joseph- 
son penetration depth Xj = \/l ± S (in normalized 
units) and the modified maximum velocity of fluxon 
""max — 1/Vl i S = c±. Sincc the fluxons in these modes 
can move with any velocity from to c±, they can not 
excite in-phase plasma waves by means of the Cherenkov 
mechanism. The state [1| — 1] with the maximum ve- 
locity c_ can not excite any plasma wave mode at all. 
Only the fluxons in the state [1|1] moving with velocity 
c- < u < c+ may, in principle, generate the out-of-phase 
mode of plasma waves. This may happen for two cou- 
pled LJJ's with different parameters. In the symmetric 
case (f) = (f)^ Cherenkov radiation does not occur since 
the solution of Eq. (O) is purely solitonic and symmetric 
while Cherenkov radiation requires the presence of the 
out-of-phase component. 

If the two LJJ's are not identical, the pure degenera- 
tion of equations does not take place and Cherenkov ra- 
diation can be generated by two fluxons moving in [1|1] 
state. In general, the derivation follows the one from 
Ref. |0| but, in Eqs. (7)-(9) of Ref. ^ the factor 2 in 
front of the sine should be changed to ( J+ 1)/2 J and the 
running coordinate to 
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III. NUMERICAL SIMULATION 

The simulations were performed using the system of 
Eqs. (M) with the simulation technique described in detail 
in Ref. |8|. Investigating Cherenkov phenomena, we took 
additional care about discretization. For simulation we 
use a finite difference method in which continuous LJJ 



is represented as a discrete chain of point-like LJJ's with 
very small spacing a <ti 1. Since the dispersion relation 
for plasma waves in a discrete medium is qualitatively 
different from the continuous one for large k, we have to 
find a parameter range in which the discreteness effects 
can be ignored. The dispersion relation for the out-of- 
phase mode of linear plasma wav£a in the coupled discrete 
array of JJ is {D' = C = J = iJB: 



>{k) = Jl + ^sin^ 
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i.e., it is periodic in k. Therefore, the fluxon dispersion 
line u! = uk intersects the dispersion relation (Hq) for any 
M > 0, i.e., even for u < c^ and, therefore, Cherenkov 
radiation can occur in simulation but not in a real contin- 
uous system. For k <^ n/a the dispersion relation (|8|) is 
identical to that of the continuous system. Since in con- 
tinuous system the radiation does not occur for fluxon 
velocities up to c_, let's calculate the wave vector kc of 
plasma wave which will be emitted in the discrete system 
at the fluxon velocity u = c_. This value is determined 
by the intersection of the fluxon dispersion line uj = c-k 
with the dispersion curve of plasma waves (fq): 



-kc = 
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Taking into account (|l^) we solve Eq. (|T^) for kc. Since 
for a <C 1 Eq. ( |lq ) is very close to the dispersion relation 
of the continuous system in the range k <^ 7r/a, it is 
sufficient to take into account small nonlinearity of sine 
in (p^), i.e. the cubic term. Solving (n9[) for kc we obtain 
the following expression: 
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Thus, for a fluxon moving with velocity m < c_, any ra- 
diation with wavelength Ac = 27r/fcc and smaller is an 
artifact of the simulation procedure. All simulation re- 
sults below were obtained with a — 0.025 and S — —0.5 
which gives Ac = 0.44 (in units of Aj). 



A. Annular junctions 

The fluxon dynamics in annular junctions was simu- 
lated using Eqs. (R) with periodic boundary conditions: 
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where N^'^ is the number of trapped fluxons in LJJ'^'^, 
respectively, and i = L/Xj is the normalized circumfer- 
ence of the stack. 

The state [1|0]. The family of I-V characteristics 
(IVC's) of the system for D' = C = 1 and different values 



of J = 0.5, 1.0, 2.0 in [1|0] state is shown in Fig. |[ The 
circumference of the stack ^ = 20 and the damping coeffi- 
cient a — 0.2 were chosen to emulate an infinite junction 
{ai ^ 1). As we predicted above, for k > 1 ( J < 1) the 
fluxon velocity exceeds c_ and Cherenkov radiation takes 
place. The phase gradient profiles for w > c_ are similar 
to those discussed below. For k < 1 ( J > 1) we find that 
Minax < c_. Notc that in this case there is a region of 
back-bending of IVC, after which LJJ^ switches to the 
resistive state. 

The presence of the back-bending on the IVC is a 
rather interesting phenomenon. Though, it is quite usual 
for non-linear resonances, in our case there is no real res- 
onance. It seems not trivial also from physical point of 
view: we pump energy into the system by means of bias 
currents but the wave (fluxon in LJJ and its image io. 
LJJ'^) as a whole moves slower. By the word "image"Ej 
we call the perturbation of the Josephson phase which 
takes place in LJJ^ due to the presence of the fluxon in 
LJJ"^. The difference between image and anti- fluxon is 
that anti-fluxon carries the flux — <I>o and image carries 
no flux. 

First, we suspected that the state of the system at the 
back-bending part of IVC is metastable and the system 
switches to some other state after a long period of time. 
However, even after waiting more than 10^ time units the 
system is still at the back bent region of the IVC and dc 
voltage on the LJJ (velocity of the fluxon) has the same 
value within the specified accuracy (better than 0.1%). 
Since in the back-bending region the state is stable, we 
suspect that there should be another reason, e.g. a struc- 
tural transformation of the solution, which may cause the 
negative differential resistance. 

We propose the following qualitative explanation based 
on numerical observations. At the point where Rd — 
du/dj approaches zero, a new strongly nonlinear solu- 
tion — a breather (fluxon-antifluxon pair) — forms in 
LJJ^. At the beginning, the distance between the fluxon 
and antifluxon in the breather is small and the amplitude 
of the breather is close to zero. The breather does not os- 
cillate since it is stretched by the bias current 7^ which 
pulls the fluxon and antifluxon in opposite directions. 
The higher the bias current, the stronger it stretches the 
breather, the larger is the distance between the fluxon 
and antifluxon. As the distance increases, the profile of 
the breather becomes similar to the profile of a fluxon and 
an antifluxon situated at some distance from each other. 
Therefore, as the distance between the fluxon and anti- 
fluxon grows from to about Aj, the friction force grows 
from to about 2Fa, where Fa is the friction force act- 
ing on the fluxon. Since the friction force grows very fast 
with 7^, much faster than the driving force, this leads 
to a back-bending region on the IVC. At some point, the 
current 7^ becomes large enough to break the fluxon- 
antifluxon pair situated in LJJ^ and, after rather com- 
plicated transient process, LJJ^ switches to the non-zero 
voltage state (R-state). 

To learn about this scenario, we simulated the tran- 



sient process of switching of LJJ^ to the R-state for 
J = 2.0. Starting from the point ^^ = 0.35 at the 
very top of the back-bending part of the IVC, the current 
was increased by a smaU amount A7 — 0.05 which was 
enough for LJJ^ to switch to the R-state. We observed a 
dissociation of the fluxon-antifluxon pair in LJJ^. Three 
snapshots of the transient process of fluxon-antifluxon 
dissociation are shown in Fig. g. 

Fig. ^ represents the IVC of the long, resonance-free 
system. The finite iength of the stack {at ^1) resuhs 
in additional effectaJ. As an example, we simulated the 
IVC of the stack with D' = C = 1, J = 0.5, S = -0.5, 
a = 0.1 and i = 7 in the [1|0] state. According to the 
results of the previous section, for this set of parame- 
ters Umax > C-, and Cherenkov radiation has to appear 
at c_ < u < Umax- The IVC is shown in Fig. 0. The 
choice of parameters was close to that of our experimen- 
tally studied sample (see below) except for the damping 
coefficient a = 0.1 which was taken larger than in the ex- 
periment in order to decrease the calculation time. For 
comparison, the IVC's for 5' = 0, and S = —0.5, J = 1 
are shown in Fig. as well. 

First, from simulations we find that for the bias current 
7 > 7* (see Fig. ^) the soliton velocity u becomes larger 
than c_, in accordance with consideration given above. 
Second, there are two steps to the right from c_. To un- 
derstand the nature of these steps we calculated the phase 
gradient profiles (f>^ ' {x) for various points of the IVC. If 
one increases 7 from to 1, the fluxon dynamics develops 
in the following way. In the region < u < c_ the fluxon 
motion is qualitatively well described by the perturbation 
approach (IS"] ^ 1)E3 and we find the fluxon in LJJ"^ and 
its image in LJJ^. Their profiles are shown in Fig. H(a) 
that corresponds to the point A on the IVC (see inset in 
Fig. H). As soon as u exceeds c_, a Cherenkov radiation 
wake arises behind the moving soliton and its image as 
shown in Fig. 0(b) which corresponds to the point B on 
the IVC. With increasing soliton velocity, the wavelength 
of the radiation A ~ 2Ti/k increases in agreement with 
(O) , and the amplitude and length of the wake quickly 
grow. The amplitude of the wake decays exponentially in 
time and space as can be seen from Fig. 0(c) correspond- 
ing to the point C on the IVC (see inset in Fig. |3|). At 
any point of the IVC the area under the soliton profile is 
J (j)^ dx = 2-K and for the image J (j)? dx = 0. 

Since the soliton moves in the annular LJJ of finite 
length £, at some velocity u the Cherenkov radiation wake 
extends over the length ~ £ so that the soliton and its im- 
age "see" their own radiation wakes after turning around 
the junction. This is illustrated in Fig. ^(d) and corre- 
sponds to point D on the IVC (see the inset of Fig. a). 
The interaction of the soliton with its Cherenkov radia- 
tion wake results in the appearance of resonances on the 
IVC of the system at u > c_ . We call these resonances 
"Cherenkov steps" (CS's). At the resonance the oscilla- 
tions of (j)"^'^ in the LJJ's induced by the moving soli- 
ton take place in phase (show constructive interference) 



with Cherenkov-generated Josephson plasma waves 0"^'^ 
[Fig. ^(d)]. Constructive interaction is possible for both 
fluxon and image at the same time since they, like out-of- 
phase plasma waves, have opposite polarity. If the fluxon 
moves in the minimum of the plasma wave in LJJ , the 
image moves in the maximum of plasma wave in LJJ^ 
and such a state is the minimum energy state for both 
the fluxon and the image. Thus, the steps CS2.3 on the 
IVC in Fig. are related to the interference between the 
moving soliton and its Cherenkov radiation. Since the 
amplitude of the emitted plasma waves decays exponen- 
tially ~ exjp{—ax), the condition of having resonances is 
a£ ^ 1. In a long system or in the system with substan- 
tial dissipation this condition is not satisfied, the inter- 
action of a fluxon with its wake does not occur and IVC 
looks smooth as shown in Fig. fl^ for J = 0.5. Our numer- 
ical study showed that for a = 0.1 individual resonances 
appear up to ^ ~ 30. Looking at Fig. n^ one can see that 
for a — 0.2 and £ = 20 resonances are not visible. 

To calculate the positions of the resonances we suppose 
that resonance occurs if there is an integer number of 
plasma wavelengths over the length of the junction, i.e., 
the wave vector k takes one of the following eigenvalues: 



2tt 
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Thus, Cherenkov radiation should lead to resonances at 
vortex velocities equal to the phase velocity (|l^) with k 
given by (pi): 
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The resonances calculated using (^ are shown in Fig. |^ 
by thin vertical lines. The density of the resonances in- 
creases with m up to infinity when u — > c_ -I- 0. Since 
the resonances have a finite width due to damping, it is 
possible to resolve only the CS's with not very large to in 
LJJ's of moderate length i. The positions of resonances 
obtained from simulation (Fig. }^) differ from the ones 
calculated using Eq. (E4) . We suggest several reasons for 
this discrepancy. First, it may result from the large am- 
plitude of the waves which become strongly nonlinear at 
the resonance. Second, the gap in the dispersion relation 
decreases because of non-zero bias current 7. To take into 
account these two factors we write down a generalization 
of Eq. d) as 



(t,'^'^{x,i) = (t)^^^ + A^'^ sin(e - Loi) 
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and substitute it into Eq. iQ) without a-term. For fucther 
calculation we will exploit the following expansionsEj: 
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and 



cos {A sin z) = Jo (A) + 2 ^ 32k (A) sin2kz 



(27) 
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We are interested only in the dc and the first ac spectral 
component of these expansions. Expanding sin((/i'^'^) as 
a sine of a sum, we use the first harmonic of ( pq ) and 
dc component of (27) to obtain the dispersion relation 
for the first harmonic of the large amplitude nonlinear 
plasma wave. This dispersion relation now depends on 
the amplitude of the plasma wave A and the bias current 
7. The resonance velocities are obtained as: 
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The average amplitude of the plasma waves A is assumed 
to be the same in both LJJ's. Expression (28) implies 
that resonant velocities Um decrease as 7 or A increase. 
This means that CS's can be bent back since both 7 and 
A increase as the bias point moves up along the CS. In 
real systems this behavior is compensated by damping 
and resonances may or may not be bent back. 

We would like to mention that the simplified expres- 
sion used in our previous workEl described positions of 
resonances rather well despite that nonlincarity of the 
waves, 7^0 and J ^ \ were not taken into account. 
This coincidence occurred because the effect of J < 1 
and the effect of the nonlincarity with 7 > have op- 
posite signs and, for our particular parameters, nearly 
cancel each other. 

The appearance of the CS's is rather convincing indica- 
tion of the Cherenkov radiation. The steps appear due to 
the constructive interaction between the fluxon and the 
image from one side and the out-of-phase plasma wave 
from the other side. The experimental observation of thft 
CS's was accomplished for the first time by our groupO 
and is described in the next section. 

The state [1 1 1] . As it was mentioned in the previous 
section, if two LJJ's are identical the [1|1] mode does not 
emit any radiation. We accomplished a series of simula- 
tions to investigate how the spread in parameters affects 
the motion of the fiuxon. Examples of numerically cal- 
culated fiuxon profiles are shown in Fig. ^. All parame- 
ters of LJJ's were identical except for the damping coef- 
ficients i? = 1.5 in Fig. |^(a) and critical currents J = 1.5 



Fig. ||(b). We found that J ^ 1 as well as R ^ \ result 
in Cherenkov emission of out-of-phase plasma waves. In- 
creasing the parameter spread one decreases the stabil- 
ity of the in-phase fiuxon step and, at some spread, the 
[1|1] state is destroyed and the in-phase fluxon step dis- 
appears. In the [1|1] state, the interaction of the fluxon 
and the image with the out-of-phase plasma wave can not 
be constructive since two fiuxons have the same polarity 
while the plasma waves have the opposite one. Thus, if 
the condition of minimum energy is held for LJJ , it is 
not held for LJJ^. As a result, the IVC has usual smooth 
behavior without additional resonances. Therefore, ex- 
perimental observation of Cherenkov radiation from flux- 
ons moving in the [1 1 1] state is a rather difficult task. 



B. Linear junctions 

The fiuxon dynamics in linear junctions was simulated 
using Eqs. (H) with boundary conditions which corre- 
spond to zero applied external magnetic field: 



,A.B 



A.B, 



C' (0-0 = C' (^,t) = o 



(31) 



In stacks of linear geometry (open ends) one may also 
observe Cherenkov radiation and resonances similar to 
that in annular stacks. In this case the fiuxon moves at 
zero applied magnetic field reflecting as antifluxon from 
the edges of the junction. In the [1|0] state, the image 
has also to change its polarity upon reflection. This mo- 
tion is similar to the fiuxon motion in single LJJ biased 
on zero-field step (ZFS). At a fiuxon velocity larger than 
c_, the Cherenkov radiation wake appears. When the 
length of the wake becomes comparable with the length 
of the junction, the Cherenkov emission forms a standing 
plasma wave. For the case of the [1|0] state the interac- 
tion between the standing plasma waves and the fiuxon is 
constructive and leads to the series of resonances (CS's) 
shown in Fig. |g. The fiuxon dynamics at the CS is rather 
similar to the fiuxon dynamics in a single LJJ at the Fiske 
Step (FS), where moving fiuxons interact with standing 
cavity mode. In the case of Fiske steps, the plasma wave 
is excited due to boundary conditions — when the fluxon 
leaves the junction, a plasma wave is emitted. In the 
case of coupled LJJ's the energy is pumped into stand- 
ing plasma wave by means of Cherenkov mechanism. The 
positions of the resonances can be calculated using the 
Eq. ( p8| ) taking 21 as the length of the junction in the 
expression for km (Pq). We have to mention that reflec- 
tion of the image results in a large energy dissipation 
that decreases the stability of the ZFS [1|0] mode (height 
of ZFS) as the coupling \S\ or the number of coupled 
junctions N increases. 

For the junctions of linear geometry one may try to de- 
tect the Cherenkov radiation directly by measuring the rf 
voltage at the junction boundary. Since the plasma wave 
oscillations are out-of-phase in the two LJJ's, it will be 



almost impossible to detect anything meaningful exper- 
imentally if the rf voltage is measured across the whole 
structure. Therefore, the sample should be designed in 
such a way that it allows to take the rf signal from one 
of the junctions. 

At the working point chosen on a CS, the frequency of 
the plasma wave is a multiple of the Josephson frequency 
^pi = TTiujj, where to is a resonance number. The spec- 
trum of the voltage oscillation (fluxon part plus plasma 
wave part) in each junction will have a higher amplitude 
of the rn-th harmonic. Nevertheless, such a behavior of 
771-th harmonic can not be considered as doubtless indi- 
cation of Cherenkov emission. Simulation shows that in 
the [1|0] state even the system with J = 1 which should 
not exhibit Cherenkov effect has enhanced amplitudes of 
high harmonics. The power is pumped into high har- 
monics because of reflection of images at the edge of the 
junction when the plasma waves are emitted in a fashion 
similar to Fiske steps. Since this plasma wave interacts 
with the fluxon, one may observe resonances on the IVC 
of the system even at m < c_ as shown by point A in 
Fig. nl We will not consider this resonance here since it 
is a subject of separate investigation. 

Thus, in the two-fold stack of linear geometry neither 
the presence of resonances on the IVC nor enhanced am- 
plitude of high harmonics in the spectra can be a reliable 
indication of the Cherenkov emission. An example of the 
spectra that correspond to the points A, B, C and D on 
the IVC (Fig. |) is shown in Fig. g. 



IV. EXPERIMENT WITH TWO-LAYER 
JOSEPHSON RING 

Experiments have been performed with (Nb-Al- 
A10a;)2-Nb annular stacks. Due to the magnetic flux 
quantization in a superconducting ring, the number of 
vortices initially trapped in an annular junction is con- 
served. The soliton dynamics can be studied here under 
periodic boundary conditions, which exclude a compli- 
cated interaction of the soliton with the junction edges. 
The mean diameter of the stacked LJ J's was D = 132 /im 
and the width W — 10 /im. The thickness of the middle 
Nb layer was about 60 nm. The normalized circumference 
of the ring at T = 4.2 K was nD/Xj = L/\j w 7, where 
\j is the Josephson penetration depth. Fluxon trapping 
was realized by applying a small bias current to the junc- 
tion while cooling the sample through the critical tem- 
— 9.2 K. Measurements were performed 



- Nb 



perature Tc 
at 4.2 K. 

In stacked annular LJJ's, trapping of a single fluxon 
is rather difficult due to the asymmetry of the required 
state [1|0]. Therefore, after every trapping attempt, the 
resulting state was carefully checked by several means. 
First, the IVC of the system was traced. A large critical 
current /c implies the [0|0] configuration. If Ic is close 
to zero, the voltage of a soliton step on IVC has to be 



proportional to the total number of vortices trapped in 
the system minus the number of pinned vortices. Second, 
the dependencies of the critical current I,, and the max- 
imum current at the single junction gap voltage Ig on 
magnetic field H applied in the plane of the sample were 
measured. The soliton configucation was deduced from 
the shape of Ic{H) and /g(i/)llJ. As shown in Fig. ||, in 
the [1|0] state Ic{H) has a minimum dX H — Q and Ig{H) 
has a maximum. The LJJ in the stack under investiga- 
tion had the ratio of critical current densities J — jf / j^ 
of about 0.5, therefore it was easy to distinguish which 
junction contains the fluxon. 

The fluxon step in the [1|0] configuration is shown in 
Fig. p. In addition to the main fluxon step, we find 
two additional steps denoted CS2 and CS3 at voltages 
Vcs2 ~ 33.3 /iV, and Vcs^ ~ 30.3 ^V. Similar IVC's 
were observed in several other stacked annular LJJ's. 

To find the experimental values for the velocities 5+ 
and c_ in the stack, we measured the FS Y^Jtage spacing 
for both in-phase and out-of-phase modesllj by applying 
a magnetic field in the plane of the sample. Using fields 
up to 35 Oe we observed two families of FS's at voltages 
above the single-junction gap voltage. In this case, one 
junction is biased at the gap voltage state and the excess 
voltage is generated by cavity resonances in the other 
jimction. We obtained Ay_ sa 30 /iV and AV+ w 56 /iV 
for the FS spacing in the out-of-phase and in-phase mode, 
respectively. From this data the coupling parameter is 
estimated as 



S 



^vl 



^vl 



AW -f- ^Vi 



-0.55 



and the limiting velocities c+ ~ 0.038c, c_ ~ 0.020c and 
Co « 0.025c are calculated. The voltage corresponding to 
the velocity cq is equal to Vb ~ 37.4 ^V. The steps CS2,3 
appear at vortex velocities above c_, and, therefore, cor- 
respond to the Cherenkov steps that are also found in 
numerical simulation (see Fig. y). The parameters for 
simulation were chosen by purpose very close to the pa- 
rameters of our sample. 

If the fiuxon is trapped in another LJJ of the stack, 
the situation J > 1 takes place and we could see the a 
small region on the top of the IVC which was bent back 
as shown in Fig. |lC|. This bending back looks very similar 
to the one found in numerical simulation and shown in 
Fig. 1^. Thus, experimental data show good agreement 
with numerical simulations for the both cases (J < 1 
and J > 1) discussed above. 



V. STACKS WITH NUMBER OF LAYERS TV > 2 

In the ideal case of high-Tj. layered superconductors the 
intrinsic LJJ's have very similar critical current densities. 
The system of equations which describes TV inductively 
coupled LJJ's was derived in Ref. and in the case of 
Resistively Shunted Junction (RSJ) model and absence of 



perturbations (without a and 7 terms, i.e. quasiparticle 
and bias currents) can be written as 



where 



ixx = S sin (f) + Scjjtt , (32) 
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is the matrix which describes the interaction between 
LJJ's, —0.5 < S" < is a dimensionless coupHng con- 
stant which for iV-fold case can be expressed as: 
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-X 



, d d 

dj smh — + 2A cosh — 

A A 



(34) 



where dj is the thickness of insulator (tunnel barrier) be- 
tween adjacent superconducting layers, d is the thickness 
of each superconducting layer and A is the London pen- 
etration depth. In the limit of high-Tc superconductors 
di,d<^ \ and 



1 



djd 



(35) 



A. Plasma Waves 

It was derived that in an TV-fold stack of LJJ's there 
are N modes of plasma waves and, accordingly, N dis- 
persion branches with A'' characteristic velocities. In the 
case of equal LJJ's parameters, the exgression for n-th 
characteristic velocity is fairly simpleEJij: 



dN) 



, in=l...N). (36) 
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Note, that for given N 



-N 



< 
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As soon as 

-N 



a fluxon is accelerated up to the velocity u > 
above the lowest Swihart velocity, Cherenkov emission of 
the slowest mode with characteristic velocity cjy may 
take place. If u exceeds several lowest Swihart veloci- 
ties, the Cherenkov radiation may consist of a mixture of 
plasma waves of different modes with different k and lu 
as illustrated in Fig. nM. An example of numerical simu- 
lation of a Cherenkov radiation wake for A^ = 7 is given 
in Ref. RtI 



B. Fluxons 

In the case of two coupled LJJ's, the asymmetry is 
essential for Cherenkov emission to appear. An impor- 
tant question which arises for the case of A^-fold stack, 
and especially for layered high-Tc materials, is whether 
Cherenkov radiation may appear in a system where all 
junctions have the same parameters. In the case A^ > 3 
there is additional decrease of the symmetry of the sys- 
tem in comparison with the case N = 2. For N = 2 and 
J = 1 both LJJ's are at the same conditions and each of 
them interacts with its neighboring junction in the same 
way. But already for TV = 3, the middle LJJ^ has two 
neighbors while LJJ's"^'*-^ have only one. Let us take 3 
coupled sine-Gordon equations and assume the [0|1|0] or 
[1|0|1] state with the symmetric solution (j)^ — cff . In 
this case one can reduce 3 equations to two equations, 
one for 6^ and another for 6^ = ch'-' : 
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(37) 

In Eq. (p^) for t/)^ there is a factor 2 in front of the cou- 
pling term. Thus, the problem of 3-fold symmetric stack 
is reduced to the problem of 2-fold asymmetric stack. 
Thus, the asymmetry appears automatically when one 
goes from A^ = 2 to iV = 3. 

The analysis presented above for the case J 7^ 1 can 
not be directly applied for asymmetric coupling term as 
in Eq. ^. We simulated Eq. (|^ for the case of [0|1|0] 
and [1|0|1] fluxon states and found that the state [0|1|0] 
is similar to the case J < 1, i.e., the IVC bends to the 

_(3) 

right in the region u > c\ which results in emission of 
Cherenkov radiation, while the state [1|0|1] is similar to 
the case J > 1 and results in the IVC with the back- 

-(3) 

bending region at m < C3 . Qualitatively, the IVC's for 
this case look like IVC's shown in Fig. |l| for J = 0.5 and 
J = 2.0. Using this result, we predict that a fluxon mov- 
ing in a layered high-Tc superconducting material causes 
Cherenkov radiation and its IVC will not have the con- 
ventional single junction asymptotic behavior. 

Since A^-fold stacks with N > 2 have lower symmetry, 
[1 1 ... 1 1] fluxon state will emit Cherenkov radiation even 
in the case of equal LJJ's parameters. If the emission is 
very strong (for A^ ^ 2), the state [1| . . . |1] can not sur- 
vive. We simulated three-fold stack in symmetric fluxon 
configuration [1 1 1 1 1] and found that this state is still sta- 
ble and emits Cherenkov radiation for u > c\ as shown 
in Fig. |l|. 

In an A^-fold stack of linear geometry, it turned out 
to be rather difficult to stabilize the motion of fluxon in 
[0|..|0|1|0|..|0] state in ZFS-like mode. Due to the pres- 
ence of images in the neighboring (relative to the LJJ's 
with fluxon) LJJ's, the reflection conditions become not 
perfect (energy dissipates in collision of images and anti- 
images at the boundary) and the fluxon tends to leave 



the junction. The ideal reflection conditions are reahzed 
in uncoupled LJJ, i.e. at 5* = 0. When \S\ and/or N 
(more images induced by one fluxon) increase, more en- 
ergy dissipates during the reflection process. For a given 
iV, at some critical value of S, the fluxon can not be re- 
flected back and leaves the system. For layered high- Tc 
superconductors (iV ^ 2) \S\ approaches its maximum 
theoretical value of 0.5. Thus, ZFS-like fluxon motion can 
most probably not be reahzed in the [0|..|0|1|0|..|0] state 
in a high- Tc layered stack with open boundaries. There- 
fore, Cherenkov radiation from a single fluxon moving in 
a stack of intrinsic LJJ's or in an A^-fold stack {N ^ 2) 
can be studied only in an artificially made annular struc- 
ture. 

In experiment, one often deals with many rather than 
with only one fluxon moving in, e.g. flux-flow mode. 
Will Cherenkov radiation appear in this case? Assum- 
ing [0| . . . |0|M|0| . . . |0] state in the hmit of dense uni- 
form fluxon chain and neglecting the junction edges it is 
possible to show that Cherenkov radiation does not take 
place since the maximum velocity of chain in each mode is 
equal to Cn . However, at high velocities when the fluxon 
chain becomes sparse, we may expect a crossover to the 
single fluxon behavior discussed above. The Cherenkov 
radiation will result in a modifled structure of the top 
part of the flux-flow branch of the IVC. 



VI. CONCLUSION 

We have shown that a single fluxon moving in a stack 
of long Josephson junctions can emit a Cherenkov radi- 
ation wake which consists of out-of-phase plasma waves. 
To generate Cherenkov radiation in the two-fold stack the 
proper asymmetry is required. In the iV-fold {N > 2) 
stack asymmetry is not necessary and a fluxon mov- 
ing faster than the lowest Swihart velocity will emit 
Cherenkov radiation. This result is of great importance 
for the physics of high- Tc materials which exhibit intrin- 
sic Josephson effect. 

For some fluxon conflgurations {e.g. [1|0]), in a system 
of moderate length {ai < 3) the constructive interaction 
of the fluxon with the Cherenkov wave results in novel 
resonances emerging on the IVC. These resonances were 
observed in experiment with annular two-fold stack and 
their positions are calculated theoretically using more rig- 
orous approach then in our previous workQ. It is shown 
that the direct observation of the Cherenkov radiation 
from the LJJ of linear geometry is problematic. 
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FIG. 1. The family of IVC's of two inductively coupled 
LJJ's for different values of J in [1|0] state. For J < 1, Umax 
exceeds c_, for J > 1, Wmax < c_. 



FIG. 2. Snapshot of the transient process of switching of FIG. 12. Cherenkov radiation of out-of-phase plasma 

the LJJ^ into R-state for J — 2.0. The current 7^ is large waves generated by fiuxons moving in the [1|1|1] state in an 
enough to break the breather. annular three-fold symmetric stack at 7 = 0.6, S = —0.5, 

£ = 20, u = 1.71 (c+ = 1.84). 

FIG. 3. Calculated dependence of velocity u on bias cur- 
rent density 7 (normalized IVC) for one soliton in a 2-junction 
ring for L/\j = 7, a = 0.1, J = 0.5, S = -0.5. The inset 
shows in detail the top part of the IVC. The IVC's for dif- 
ferent configurations are shown for comparison: dotted line 
corresponds to single annular junction (5* = 0), thin line to 
the stack in the [1| — 1] state. 

FIG. 4. Profiles of 0j'^(x) in both LJJ's. (a), (b), (c), 
(d) correspond, respectively, to IVC points A, B, C, D shown 
in the inset of Fig. 

FIG. 5. Cherenkov radiation of out-of-phase plasma waves 
generated by two fluxons moving in [1|1] state in annular 
two-fold stack at 7 = 0.6, S = -0.5, £ = 20. (a) R = 1.5, 
u = 1.27; (b) J = 1.5, u = 1.34. 

FIG. 6. The first ZFS with Cherenkov steps on the top of 
it on IVC of linear two-fold stack. 

FIG. 7. Spectra of rf voltage on LJJ , on LJJ^ and on 
both junction (black, white and gray) for different bias points 
are shown in plots (a)-(c), respectively, and correspond to the 
bias points A (7 = 0.16), B (7 = 0.2}, C (7 = 0.35) and D 
(7 = 0.53) on the IVC shown in Fig. 

FIG. 8. Dependence of the critical current 7c and the 
maximum current at the gap voltage Ig of the magnetic field 
H. The shape of these characteristics implies [1|0] fluxon 
state. 

FIG. 9. Experimental IVC of the double-layer annular 
junction in the [1|0] soliton configuration for J « 0.5. The 
Cherenkov resonances are marked CS2,3. The inset shows the 
geometry of the sample. 

FIG. 10. IVC of the double-layer annular junction in the 
[1|0] soliton configuration for J « 2.0. The back-bending 
region is visible. 

FIG. 11. Dispersion relation of 7-fold stack with coupling 
constant S = —0.5 and dispersion line of a fluxon. At velocity 
u > c„_3 3 modes of plasma waves are emitted. 
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